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Abstract 

We study the tick dynamical behavior of the bond futures in Korean Futures Exchange (KOFEX) 
market. Since the survival probability in the continuous-time random walk theory is applied 
to the bond futures transaction, the form of the decay function in our bond futures model is 
discussed from two kinds of Korean Treasury Bond(KTB) transacted recently in KOFEX. 
The decay distributions for survival probability are particularly displayed stretched expo- 
nential forms with novel scaling exponents (3 = 0.82(KTB 203) and (3 = 0.90(KTB112), 
respectively, for our small time intervals. We obtain the scaling exponents for survival prob- 
ability e = 17 and 18 decayed rapidly in large time limit, and our results are compared with 
recent numerical calculations. 
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1. INTRODUCTION 



Recently, the investigation of economic phenomena between economists and physicists 
has received considerable attention as one interdisciplinary field. 1 ' 2 Among many outstanding 
topics, there has been mainly concentrated on price changes in open market, 3 ' 4 the distribu- 
tion of income of companies, the scaling relation of company's size fluctuations, the financial 
analysis of foreign exchange rates, 5 " 7 and tick data analysis of bond futures. 8-10 Particu- 
larly, the study of fluctuations in economics can lead to a better understanding for the scaling 
properties based on methods and approaches of physics. The power law distribution, the 
stretched exponential distribution, and Levy-like distribution have been really elucidated 
the general mechanism in diverse fields of science and technology. Mantegna and Stanley 4 
reported that price changes in a stock market scales as a power law of the index change 
with a scaling exponent a = 1.4. Takayasu et a/. 5 have showed the stretched exponential 
distribution from results of analysis observed by tick data of yen-dollar exchange rates and 
have recently investigated extensively several hyper-inflations 11 in which the rapid growth of 
price changes is approximated by a double exponential function. Lee at al. 12 have obtained 
the fluctuation in the gross domestic product that the width of the distribution scales as a 
power law of the gross domestic product with a scaling exponent (5 = 0.15. Otherwise, many 
works 5 ' 6 ' 7,13,14 have reported that the price fluctuations follow the anomalous power law from 
the stochastic time evolution equation, i.e. Langevin-type equation. 

In the stochastic process, one of the well-known popular subjects for the regular and dis- 
ordered systems is the random walk theory, 15 and this theory has extensively developed to 
the continuous-time random walk theory which is essentially described both by the transition 
probability dependent of the length between steps and by the distribution of the pausing 
times. 16 Until now the transport phenomena for the motion of tracers have largely been ex- 
tended to the reaction kinetics, 17,18 the strange kinetics, 19 fractional diffusion equations, 20,21 
and generalized Langevin and Fokker-Planck equations. 22-26 
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Recently, Scalas 8 presented the correlation function for bond walks from the time series 
of BUND and BTP(Buoni del tesoro Poliennali) futures exchanged at the London Interna- 
tional Financial Futures and options Exchange(LIFFE). Scalas et al. 9 have discussed that 
the continuous-time random walk theory, formerly introduced by Montroll and Weiss, 16 , is 
applied to the dynamical behavior by tick-by-tick data in financial markets. Mainardi et 
al. 10 have also observed the waiting-time distribution for bond futures traded at LIFFE, 
London. In this paper, theoretical and numerical arguments for the volume of bond futures 
traded at KOFEX of Pusan are presented. Mainly, we consider the number of transactions 
for two different delivery dates and find the decay functions for survival probabilities in our 
bond futures model. Sec. 2 is devoted to the survival probability in the continuous-time 
random walk theory and to the decay function of survival probabilities. In Sec. 3, we reports 
numerical results and some concluding remarks. 



2. THE SURVIVAL PROBABILITY AND THE STRETCHED EXPONEN- 
TIAL FUNCTION 



First of all, we will focus on the continuous-time random walk theory 16 in order to discuss 
the survival probability. Since the pausing-time distribution ip(t) between two consecutive 
transactions and the price l(t) can be modelled as random variables in financial markets, we 
consider that the coupled transition probability density ^f(l,t) for jump I and pausing-time 
t is defined by 

V(l,t)=p(l)iP(t) (1) 

where p(l) is the transition probability dependent of the length between steps and ip(t) is the 
pausing-time distribution. Summing over all / with the periodic boundary condition, Eq.(l) 
is given by 

E*(M) = EpW(0=iW)- (2) 

I I 

Since R n {l, t) is the probability density arriving immediately at lattice point I at time t after 
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one random walker goes to n steps, R n (l,t) is satisfied with the recursion relation as 

Rn(l,t) = ]T f dt'p(l-l')i;(t-t')R n ^(r,t'). (3) 

In order to find the probability density, P(l,t), which is existed at lattice point I at time t, 
we introduce the relation between P(l,t) and R(l,t — t') such as 

P(l, t) = f dt'4>(t')R(l, t-t'), (4) 
J o 

where 

oo 

R{l,t) = Y,Mt)P(l,n)- (5) 

n=0 

The survival probability, ^{t), which is stayed for the time t after arrivng at an arbitrary 
lattice point, can be expressed in terms of 

V(t) = 1 - f dt' V^(M') = 1 - f ' dt'W). (6) 
J o \ Jo 

Hence, after substituting Eq. (6) into Eq. (5), the Fourier-Laplace transform of the proba- 
bility density, i.e. Montroll- Weiss equation, is described 16 as 

P(q,u) = Y dtP(l,t)exp(-iql-ut) 
R(q,u)[l-^(u)] 



u 

P(q,0) 



(7) 



[1 -*(<?,«)] u 

where the generalized structure function, i.e. the Fourier-Laplace transform of Eq. (1) is 
given by 

V{q,u) =p{q)^j{u) (8) 

and P(q, 0) is the Fourier transform of the initial condition P(l, 0). 

Let us introduce the following scaling behavior in the long jump and long time limits: 

P(q) ~ 1- I Q \ a ,Q^ 0,0 < a < 2, (9) 

and 

ip(u) ~ 1 -u^,u-> 0,0 < < 1. (10) 



We can also present approximately as 

p(q) = exp(- | q | Q ),0 < a < 2 (11) 



and 



^(u) = t ^,0<(3<1, (12) 



where Eq. (11) means one type of the characteristic function for Levy stable distribution. 
Hence the inverse Fourier transform of Eq. (11) is found to scale approximately with p(l) ~ 
| / for < a < 2 in the large / limit. 

Since the Laplace transform of survival probability from Eq. (6) is given by 

m = iz|M, (13) 

LI 



we find that 



u p-i 



after substituting Eq. (12) into Eq. (13). The survival probability, i.e. the inverse Laplace 
tranform of Eq. (14), can be derived as 

= E P (-1*), (15) 
where the Mittag-Leffler function of order ft 27 ' 28 is 

oo j./3n 

E ^=U- lT m^ry (16) 

We also note that the decay distribution of the survival probability appears in the fractional 
relaxation process for < (3 < 1 while it includes to the fractional oscillation process in the 
case of 1 < P < 2. 

According to Eq. (16), the survival probability ty(t) for sufficiently small times reduces 
to a stretched exponential: 

{ltf n 



E^tf] ~ 1- 



T((3n+1) 



^ exp| "WTi)l- (17) 



where 7 is a time-scale factor. In the large t limit, the survival probability fy(t) from remained 
transaction quantities less than 10 percentages of total tick data scales as a power law 

#(f) elogt, (18) 

where e is the scaling exponent for large times. In our scheme we will make use of Eq. (16)- 
(18) to find out the dynamical behavior for the survival probability, and these mathematical 
functions would be expected to lead us to more results. 

3. RESULTS AND CONCUSIONS 

In order to analyze survival probabilities for bond futures with delivery date, our un- 
derlying asset is a Korean Government bond traded at KOFEX. Since the continuous-time 
random walk theory can sucessfully account for the time evolution of financial markets that is 
non-local and non-Markovian, we concentrate on the survival probabilities ty(t) independent 
of transaction places. We here consider two different delivery dates : December KTB112 
and March KTB203. The tick data for December KTB112 were taken from July 2001 to 
December 2001 while we used the tick data of March KTB203 transacted for six months 
from October 2001. The KOFEX opens for 6 hours per one day except for weekends, and 
the total tick data of transactions are, respectively, about 2.3 x 10 5 (KTB112) and 2.1 x 10 5 
(KTB203). 

Figs. 1 and 2 depict, respectively, the survival probability ty(t) as a function of time step 
t for December KTB112 and March KTB203 delivery dates. We will assume that one time 
step is the transaction time evoluted for 10 seconds, and total time steps are approximately 
2.9 x 10 5 (KTB112) and 2.5 x 10 5 (KTB203), neglecting time intervals of the market closure. 
The slope of the solid line from Fig. 1 gives approximately the streched exponential function 
of Eq. (17) with j3 = 0.90 and 7 = 1.5 x 10~ 5 for December KTB112 delivery date while 
the survival probability ^(t) decays with (3 = 0.82 and 7 = 2.0 x 10~ 5 from Fig. 2 of March 
KTB203 delivery date. Figs. 3 and 4 is also the plots of log ty(t) as a function of log t in Eq. 
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(18), by taking continuous tick data as shown in Figs. 1 and 2. On the other hand, we can, 
respectively, obtain approximately the scaling exponents e = 17(KTB112) and 18(KTB203) 
decayed rapidly in large time t limit, as plotted in the dashed line of Figs. 3 and 4. 

In conclusion, we have presented the dynamical behavior of the continuous tick data 
for our bond futures, based on a continuous-time random walk theory recently proposed by 
Scalas et a/. 9,10 For two kinds of Korean Treasury Bond(KTB203 and KTB112) transacted 
recently in KOFEX, the decay distribution for survival probability in our small time intervals 
is displayed stretched exponential form while it scales as a power law for large times. 

In Scalas model, they introduce BUND futures with two different delivery dates : June 
1997 and September 1997. The survival probabilities for their BUND futures have the form 
of the well-known Mittag-Leffler function interpolated between two limitting behaviors: the 
stretched exponential function exp[— (yt) /T(f3 + 1)] for small time intervals, and the power- 
law function (^t) -13 /Y{1 — (3) for larger time intervals. The survival probability has an index 
P = 0.96 and a time scaling factor 7 = 1/12 for BUND futures of June 1997, in Figs.l — 4 of 
Ref. [10], while it has (3 = 0.95 and 7 = 1/12 for BUND futures of September 1997. As shown 
in Figs. 1 and 2, we have the index (3 = 0.90 and a time scaling factor 7 = 1.5 x 10~ 5 for 
KTB112 and j3 = 0.82 and 7 = 2.0 x 10" 5 for KTB203 in our small time intervals. However, 
our key results for the index (3 arenot in good agreement with those of Scalas model. It is 
in practice found that our scaling exponents (3 are smaller than those of Scalas model 8,9 , 
because their transaction quantities are slightly less than ours. As transaction quantities 
increase more and more, we would suggest that the scaling exponent (3 take a value near 
\. Hence, when (3 is varied in the range \<f3<f3 c = 1, it can be available for several future 
models to carry out the continuous-time random walk. Here (3 C is the critical exponent 
which has the so-called marginal value of (3 from the dynamical behavior obtained in Scalas 
model. If transaction quantities decrease less than those of Scalas model, then the value (3 
becomes a larger value than one, and for this case the curve of the survival probability may 
be decreased rapidly. 
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In future, we will present the correlations for bond walk from the tick data of bond 
futures in Korean financial markets and in detail compare our results with bond futures 
transacted in other nations. It is also expected that the detail description of the continuous- 
time random walk theory will be used to study extensions of foreign financial analysis of 
won(Korea)-yen(Japan) and won(Korea)-dollar(USA) exchange rates. 
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FIGURE CAPTIONS 



Figure 1. Plot of the survival probability for December KTB112 as a function of the time t. 
The solid line is the stretched exponential function for (3 = 0.90 and 7 = 1.5 x 10~ 5 . The 
continuous tick data for December KTB112 were taken from July 2001 to December 2001. 

Figure 2. Plot of the survival probability for March KTB203 as a function of the time t. 
The solid line is the stretched exponential function for (3 = 0.82 and 7 = 2.0 x 10~ 5 . The 
continuous tick data for March KTB203 were taken from October 2001 to March 2002. 

Figure 3. Log-log plot of the survival probability of the time t for December KTB112. The 
solid line is the stretched exponential function for f3 = 0.90 and 7 = 1.5 x 10~ 5 , and the 
dashed line is scaled approximately as a power law with a scaling exponent e = 17. 

Figure 4. Log-log plot of the survival probability of the time t for March KTB203. The solid 
line is the stretched exponential function for (3 = 0.82 and 7 = 2.0 x 10~ 5 , and the dashed 
line is scaled approximately as a power law with a scaling exponent e = 18. 
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